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Abstract 

We find Marginal Fermi Liquid (MFL) like behavior in the Hubbard model on 
a square lattice for a range of hole doping and on-site interaction parameter 
U. Thereby we use a self-consistent projection operator method. It enables us 
to compute the momentum and frequency dependence of the single-particle 
excitations with high resolution. The Fermi surface is found to be hole-like in 
the underdoped and electron-like in the overdoped regime. When a compar- 
ison is possible we find consistency with finite temperature quantum Monte 
Carlo results. We also find a discontinuous change with doping concentra- 
tion from a MFL to Fermi liquid behavior resulting from a collapse of the 
lower Hubbard band. This renders Luttinger's theorem inapplicable in the 
underdoped regime. 

PACS: 71.10.-w, 79.60.Bm, 71.18.+y, 74.72.-h 
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After the discovery of high-temperature superconductivity in layered Cu based per- 
ovskites it was found that those materials exhibit quite unusual properties in the normal 
state. For example, in underdoped materials, i.e., for hole concentrations less than the one 
leading to the highest superconducting transition temperature, the temperature dependent 
resistivity is found to be p(T) ~ T in the normal state. Also the nuclear relaxation rate, 
e.g., of YBa 2 Cu 3 07 has an unusual temperature independent contribution. 

Aiming for an explanation of these strong deviations from a normal metal behavior Varma 
et al. 1 developed the Marginal Fermi Liquid (MFL) theory. This theory assumes that the 
frequency 00 and temperature T dependent self-energy S(cj,T) of the electrons behaves for 
uj > T like ReS(w,T) ~ cjln|cj| and ImS(cj,T) ~ in contrast to ordinary Fermi liquid 
theory where ReS(o;, T) ~ uj and Im£(u;, T) ~ uj 2 holds. Note that at zero temperature the 
MFL form of the self-energy implies a diverging effective mass at the Fermi energy. With 
these assumptions most of the observed strong deviations from normal metal behavior could 
be explained surprisingly well. However, the microscopic origin of MFL behavior of the 
self-energy has remained an open problem. 

There have been detailed studies of the two-dimensional (2D) Hubbard model as a sim- 
ple model for the high-T c cuprates 2 mainly by using advanced numerical techniques. We 
mention in particular the Lanczos method 2 , the Quantum Monte-Carlo (QMC) 3 ~ 5 method 
or calculations based on the Dynamical Cluster Approximation (DCA) 6 . 

A perturbation analysis of the half-filled case at T = has shown that in the weak 
Coulomb interaction limit a MFL type of self-energy is obtained 7 . It is due to the van 
Hove singularities which one is dealing with in this particular case. But we know that 
electron correlations are strong in the superconducting cuprates and that at T = the 
system is an antiferro magnet 2,5 ' 8 . It is also known that by hole doping the antiferromagnetic 
correlations are rapidly suppressed 2 ' 5 ' 6,9 . Nevertheless MFL behavior continues to exist in 
the underdoped regime and the question is whether or not it can be explained within the 
2D Hubbard model with fairly strong interactions. 

The aim of this letter is to demonstrate that MFL behavior can indeed be derived from 
a doped 2D Hubbard model on a square lattice at T = and large on-site interaction. This 
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has become possible with the help of a recently developed self-consistent projection operator 
method (SCPM) 10 . It allows for high resolution calculations of the self-energy as regards 
its momentum and energy and avoids certain problems previous numerical calculations have 
had to face. 

The SCPM is an extension to the nonlocal case of a projection operator coherent poten- 
tial approximation 11 . The latter was shown to be equivalent to the many-body CPA, the 
dynamical CPA as well as the dynamical mean- field theory 12 . In the following we outline 
briefly the main equations which are used before we present the numerical results demon- 
strating MFL behavior. More detailed derivations of the equations are found in the original 
literature 10-12 . 

Starting point is the retarded Green function 

G k( z ) = 1 a 7 n • (!) 

Here z = uj + iS where 5 is a positive infinitesimal number, is the Hartree-Fock one- 
electron dispersion measured from the Fermi energy and (z) is the self-energy calculated 
from the nonlocal memory matrix according to 

k k (z)=U 2 J2M j0 (z)eMik-Rj) ■ (2) 

3 

While U denotes the Hubbard on-site interaction, Rj is the position vector of site j. We 
calculate M^z) by using an incremental cluster expansion up to two-sites contributions 

M u (z) = Mjp(z) + J2 (MP(z) - M^{z)) , (3) 



M^(z) = M^(z), (4) 

where M^\z) and M^(z) are matrix elements of the cluster memory matrices M^(z) 

- ( c ) 

(c = The latter can be expressed in terms of a "screened memory matrix" M (z) 

and a matrix lS°\z) which describes on-site excitations. It is 



M (c) (\ - L<«> • M (C)N 



(5) 

lm 
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The matrices have dimensions lxl when c = % and 2x2 when c = (ij). Specifically 
L^(z) = C/(l — 2(n i _ cr ))/[(nj_ cr )(l — (rij-o-))] while (z) is a diagonal matrix with elements 



screened memory matrix is calculated from renormalized perturbation theory as 



with An = [<n i _ -)(l - (n i _ ff ))]/[(n i _ (7 ) c (l - (tw)^] and A^- = 1. Here (n ia ) c = 
J deplf (e) / '(e) with /(e) denoting Fermi's distribution. The matrix pfj{e) describes the 
density of states of a system with an empty site i (or sites i and j) embedded in a 
medium with a coherent potential Tl(z). This coherent potential is determined self- 
consistently from E(z) = N^^-^k^k^ wnere N * s the number of sites. Moreover, 



X(e,e',e") = /(-e)/(-e')/(e") + /(e)/(e , )/(-e"). We want to emphasize that in Eqs. (2)-(4) 



all memory matrices with site % separated sufficiently far from site j are taken into account 
until convergency is obtained. 

It follows from Eqs. (5) and (6) that the self-energy reduces in the limit of small U 
to second order perturbation theory while in the limit of large U the exact result of the 
atomic limit is reproduced 10 . Although the above computational scheme looks at first sight 
somewhat difficult to handle, this is not really the case. In fact, it allows us to calculate the 
self-energy directly without having to do a numerical analytic continuation or an interpola- 
tion in k space. Therefore we obtain for it a high resolution in energy and momentum. In 
the numerical calculations we have done, we assumed a paramagnetic ground state since the 
antiferromagnetism is suppressed away from half-filling. Whenever possible we have made 
comparisons with Quantum Monte-Carlo (QMC) results and the agreement was always very 
satisfactory. 

In Fig. 1 the momentum dependent excitation spectrum is shown in the underdoped 
case for U = 8 (in units of the nearest-neighbor transfer integral) and T = 0. One notices 
an empty upper Hubbard band centered around the M point and a flat quasiparticle band 
crossing the Fermi energy ep. There is also incoherent spectral density near the T point re- 
sulting from the lower Hubbard band. Also shown are QMC results for finite temperatures 5 . 
Results for the Fermi surface are shown in Fig. 2. For n = 0.95 (underdoped case) a hole- 



L^\z) and lP\z). As usual (n i<T ) is the average electron number at site % with spin a. The 
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like Fermi surface is obtained. Due to a collapse of the lower Hubbard band the portion of 
the flat band around the X points sinks below the Fermi level. In the overdoped regime 
the Fermi surface is electron-like. It is seen that Luttinger's theorem 13 does not apply here. 
These results agree with the ones obtained from the DC A 6 . 

By taking numerical derivatives of A^, (z) we have determined the momentum dependent 
effective mass m^ = 1 — dReAfc(0 + ) / du in the underdoped regime (see Fig. 3). To our 
knowledge this is something which could not be done before. For doping less than 2 % my, 
changes strongly between the minimum value at the M point and the maximum value at 
the X point, while for dopings larger than 2 % the momentum dependence of mp, is weak 
with a maximum at V and not at X as in the underdoped regime. Most important is the 
strong dependence of my, near the X point on the chosen step size 5oj when the derivative 
is taken. This is a clear signature of MFL behavior. Because of the numerical derivative 
taken of ReA^ at uj = we obtained a finite value of my, ~ \n5u instead of a divergency. In 
fact, in the limit of vanishing hole doping (half-filled case) we find a typical MFL behavior 
of Aju(z) for U — 8 and k = (ir/2,ir/2) like in the weak interaction limit. This is shown in 
Fig. 4 where the two cases are compared. Because of these features we conclude that for 
U — 8 and doping less than 2 % MFL theory applies, while for doping concentrations of 
more than 2 % normal Fermi liquid theory is valid. 

The different nature of the states for doping concentrations 5h < 0.02 and 5h > 0.02 is 
clearly seen in the density of states (DOS) presented in Fig. 5 for U — 8. Consider first the 
case 5 h < 0.02. With increasing doping concentration spectral density is shifted from the 
lower to the upper Hubbard band, or more generally to higher energies. As a consequence 
the peak in the DOS remains at ep, i.e., it does not shift for small doping concentrations. 
Therefore the self-energy has to good approximations the same frequency dependence as 
for half-filling. This is the origin of MFL behavior. When 5h > 0.02 the lower Hubbard 
band has essentially collapsed and the peak in the DOS moves away from ep. In that case 
conventional Fermi liquid behavior sets in. This occurs in a rather discontinuous way. This 
is seen from the phase diagram shown in Fig. 6. The MFL regime with a still existing lower 
Hubbard band is separated by a region in which two self-consistent solutions are found from 
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the one with a collapsed lower Hubbard band. One may consider that transition as one from 
quasi-localized electrons to fully itinerant ones. Note that when 5 h < 0.01 the MFL regime 
is divided into two regimes. The two solutions found are only slightly different in the amount 
of the reduction of the lower Hubbard band. For fixed value of Sh the discontinuous behavior 
of the self-energy as function of U implies that Luttinger's theorem 13 is not applicable. For 
U < 6.5 the MFL behavior at half-filling changes smoothly to a Fermi liquid state. 

In summary, by using the SCPM we could calculate the zero-temperature self-energy 
h-fe(z) for the 2D Hubbard model with high resolution with respect to uo and k. We find a 
MFL like behavior for quite a large range of hole doping concentration and U. We obtain 
there a strong momentum dependence of the effective mass. In cases where a comparison 
with finite temperature results 3-6 can be made the agreement is good. When U > 6.5 a 
discontinuous change takes place with increasing hole concentration from more localized 
electrons with a lower Hubbard band to fully itinerant ones. In the latter case the lower 
Hubbard band is absent. It is precisely the transfer of spectral density from the lower 
Hubbard band to higher energies which results in MFL behavior at low hole concentrations. 
Above the upper discontinuity lines Luttinger's theorem is not applicable. Very close to 
half-filling long-ranged antiferromagnetic correlations are expected to modify the present 
results. Those correlations have been neglected here. 
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FIGURES 

FIG. 1. Single-particle excitation spectra along high symmetry line for U = 8 and T = in 
unit of the nearest-neighbor transfer integral: electron occupation number n = 0.95. Open circles 
with error bars are the QMC results 5 at T = 0.33. The dashed curves show the Hartree-Fock 
contribution e&. 

FIG. 2. Excitation spectra at the Fermi energy showing Fermi surfaces for n=0.95 and 0.80. A 
80 x 80 mesh was used for calculations in the Brillouin zone 

FIG. 3. Momentum-dependent effective mass as a function of doping concentration. Closed 
circles: maximum value of at X(ir,0) for the hole concentration 5h = 1 — n < 0.02 and T (0,0) 
for 0.02 < 5h, open circles: average m^, closed triangles: minimum value at M(tt,tt). Numerical 
derivatives are taken with respect to energy fraction 5u> = 0.05. For the maximum mk, also 
results for 5uj = 0.005 are shown (+). The momentum-independent effective mass in the single-site 
approximation (SSA) is shown by the dashed curve. 

FIG. 4. Real part (thin solid line) and imaginary part (solid line) of the self-energy at the k 
point (7r/2, 7r/2) and half filling. Corresponding results of second-order perturbation theory are 
shown by dashed lines. 

FIG. 5. Densities of states (DOS) for the two self-consistent solutions at n = 0.98. One 
n = 0.98+ (solid curve) is smoothly connected with the region n > 0.98, while the other n = 0.98— 
(dot-dashed curve) with the region n < 0.98. The DOS for n = 1 is also shown (dotted line). Note 
that 5 = 0.02 (imaginary part of energy) was used in the numerical calculations, so that the peak 
at uj = remains finite for n = 1. 

FIG. 6. Phase diagram showing discontinuity lines. The different transition lines are explained 
in the text. The regimes with two self-consistent solutions are shown by hatched areas. Dashed 
lines are extrapolations. 
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